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ABSTRACT 

Context. Even though it is broadly accepted that single Be stars are rapidly rotating stars surrounded by a flat rotating circumstellar 
disk, there is still a debate about how fast these stars rotate and also about the mechanisms involved in the angular-momentum and 
mass input in the disk. 

Aims. We study the properties of stars that rotate near their critical-rotation rate and investigate the properties of the disks formed by 
equatorial mass ejections. 

Methods. We used the most recent Geneva stellar evolutionary tracks for rapidly rotating stars that reach the critical limit and used a 
simple model for the disk structure. 

Results. We obtain that for a 9 Mq star at solar metallicity, the minimum average velocity during the Main Sequence phase to reach 
the critical velocity is around 330 km s _1 , whereas it would be 390 km s _I at the metallicity of the Small Magellanic Cloud (SMC). 
Red giants or supergiants originating from very rapid rotators rotate six times faster and show N/C ratios three times higher than those 
originating from slowly rotating stars. This difference becomes stronger at lower metallicity. It might therefore be very interesting 
to study the red giants in clusters that show a large number of Be stars on the MS band. On the basis of our single-star models, we 
show that the observed Be-star fraction with cluster age is compatible with the existence of a temperature-dependent lower limit in 
the velocity rate required for a star to become a Be star. The mass, extension, and diffusion time of the disks produced when the star 
is losing mass at the critical velocity, obtained from simple parametrized expressions, are estimated to be between 9.4 • 10~ 12 and 
1.4 • 10~ 7 M Q (3 • 10~ 6 to 4.7 ■ 10~ 2 times the mass of the Earth), 2000 and 6500 R Q , and 10 and 30 yr. These values are not too far 
from those estimated for disks around Be-type stars. At a given metallicity, the mass and the extension of the disk increase with the 
initial mass and with age on the MS phase. Denser disks are expected in low-metallicity regions. 

Key words, stars: general - stars: evolution - stars: rotation - stars: Be - stars: mass-loss 



J> | 1. Introduction From an obse r vational point of v i ew, some authors 

k> ^, • , , • , „ ... (iFremat et all 12005b ICranmerl (20051 120091 iHuang et al.l 1201 Oh 

^ , There is general consensus that classical Be stars are mtnnsi- found ^ fie ^ in rticular ear i y . type Be stars, seem to 

b ; cal y rapidly rotating B-type stars, with analogs among O and fee rQtati bdow their critical ^ ^ fa ^ SQme cases ^ 

■ ■ early- A types (a lthough they are much less frequent, see e.g. determination of rotationa i velocities could be underestimated 

| Zorecetal.| | 2005|) , surrounded by a gaseous decretion viscous (Xownsend et al 2004) Also> Rivinius et al (2001) found for 

disk. Such a disk, geometrically thin and in nearly Keplenan some earl t Be stars such as Cen that the beati of 

rotation, can explain many of the c haracteristics that define Be multi iodic non . radia i pu l sa tions seems to be linked with the 

stars (see for e xample the reviews by| Porter & Rivinius | 2003 |and outbu rsts. Nevertheless, the lack of observational evidence 
Carciofi 2011). 



for large-amplitude pulsations or other strong mechanisms that 

However, there are still many open questions regarding Be couW hd tQ driye lar momentum and mass int0 a circum . 

stars: do they reach the critical limit, or do they rotate at 70-80% steUar difjk in many of ±e ^ fayours the idea of weaker 

of their critical limit? Does binanty play a role in the appearance ses that would be able to overcome the effective gravity 

of the Be phenomenon? Which are the mechanisms involved in and M ±e matter g from ±e hoto here int0 the disk: 

the angular-momentumand mass input in the disk? As noted thus rotation dose tQ the critical Hmit ifj expected ( Howarth 

byl Stee&Meilland l (|2009|), it seems that even if we consider 2QQ7) This idea is also s rted b interferometric 

the single stars that present the Be phenomenon, they do not all observationSi which have allowed to determine the oblateness 

share a common mass-ejection mechanism. 
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of some Be stars (Domiciano de Souza et al] 20031: Zhao et al 



2011 



iTouhami et alJ 1201 it iMeilland et al.l 120121: iKraus et al. 



2012), suggesting that some of them are close to their critical 
rotation. 

Interestingly, theoretical studies on stellar evolution show 
that along the main-sequence (MS) evolution, the single 
intermediate-mass stars with a sufficiently high rotational veloc- 
ity on the zero-age main sequence (ZAMS) can reach equatorial 
velocities near the critical valu^J thanks to the transfer of an- 
gular mo mentum from the inner contracting par ts to the outer 
reg ions dMevnet et alj|2007t lEkstrom etai1l2008l) . 

lEkstrom et al. (2008) studied how stars with different 
masses, different metallicities, and different initial rotational ve- 
locities evolve towards the critical limit. In that work, the com- 
putation stopped when the model reached the critical velocity. In 
the study we present here, the following c hanges and improve- 
ments with respect to Eks trom et al.l (120081) were made: 

- The numerical method for computing the surface velocity 
as resulting from the various transport proces ses and mass 
loss w as improved (as explained in detail in Georg v"et al.l 
20121 hereafte r Pape r I). In contrast to the computation by 



Ekstro m et al. (2008), in which the deformation of the star 



and the angular momentum contained in the envelop^ were 
neglected, we account here for both effects. The envelope 
contains one ten-thousandth of the total stellar mass during 
the MS phase. We assumed that this region has the same an- 
gular velocity as the first layer of the interioiQ. This devel- 
opment allows for a much better tracing of the evolution of 
the total angular-momentum content of the star. As a result, 
our models now reac h the critical limit less easily than in 
lEkstrom et al.l d2008l) . as discussed in Sect. [3] In addition, 
the effects of anisotropi c winds were accou nted for in the 
computations following Georgv et all (1201 ll) . 

- The computations were now pursued all along the MS phase 
until the helium flash for stars with M m \ < 2 Mq, the early 
asymptotic giant branch for 2.5 Mq < M m \ < 9 Mq or the 
end of core-C burning when M m \ > 9 Mq, even for the mod- 
els that reach the critical limit. This allows us to explore the 
consequences of fast rotation beyond the MS phase. 

- In the present w ork, the physical ing redients are the same as 
in the paper by Ekstrom et al.l (120121) . and differ from those 
used by Ekstr om et al.l (12008). as discussed in Sect. [3] 

Using this new version of the Geneva stellar evolution code, 
our group has produced grids of stellar-evolution tracks at solar 
metallicity (Z = 0.014) in the mass range of 0.8 to 120 Mq, both 
rotating and non-rotating (Ekstrom et al. 2012), and showed that 
the models rotating with Vini/Vcrit =0.4 provide a reasonable fit 
to important observed characteristics of stars. In Paper I we pre- 
sented a grid to explore in detail the mass domain between 1 .7 
to 15 Mq, which corresponds to early- A, -B, and late-O spectral 
types, at different metallicities and with initial angular velocity 
between < £2j n j/£2 CI i t < 0.95. In this last paper, the effects of 
rotation on the tracks in the HR diagram, the lifetimes, the sur- 
face velocities, and abundances were discussed. 

Here we specifically focus on a sub-set of these models: 
those whose surface velocity at the equator reaches the critical 

1 By critical velocity we mean here the velocity required for the cen- 
trifug al force to counterbalan ce the gravity at the equator of the star (see 
e.g. lMaeder & Mevnetll200Ch . 

2 The envelope is the region where convection is treated in a non- 
adiabatic way and ionisation may not be complete. 

3 The stellar interior is defined as the region where convection is adi- 
abatic and ionisation is complete. 



limit during the MS phase. With the study of these models, we 
aim to address the following questions: what can be said about 
the mechanical mass-loss rates when the star reaches the critical 
limit? What would be the extension and lifetime of a stationary 
circumstellar viscous decretion disk formed by matter expelled 
from the star? What is the mass-loss rate expected from the star- 
disk system? How do the properties of the disk vary as a function 
of the mass, metallicity and age of the star? To find the an swer, 
we used the analytical developme nt in IKrticka et at] d201 lb to- 
gether with the results presented in Pap er ll 

The paper is organised as follows: the main physical ingredi- 
ents of the stellar models are briefly recalled in Sect. |2] In Sect. [3] 
we discuss the properties of models that reach the critical veloc- 
ity during their MS phase: average values for the rotation ve- 
locity, properties of red giants and supergiants originating from 
very rapid rotators, and the frequency of rapid rotators in clus- 
ters of various ages. We discuss the mechanical mass-loss rates 
obtained when the models reach the critical limit in Sect. |4] 
Section [5] discusses the properties of the disks around models 
that rotate at t he cri tical limit using the disk model described in 
IKrticka et alJ yOl 1). The possible relations between stars at the 
critical limit and classical Be-type stars are discussed in Sect. [6] 
Finally, the main conclusions are presented in Sect. Q 

2. Stellar models 

The stellar models analysed here have masses between 1 .7 and 
15 Mq, and metallicities of Z = 0.002 (Z SM c), 0.006 (Z, mc), and 
0.014 (Zq). They are a sub-set of the grid presented in iPaper fl 
since we consider here only the models that reach the critical 
limit during the MS phase, i.e. models with a> m i = Q; n i/Q cr i t > 
0.80. A detailed de scription of the physics included in the mod- 
els can be found in Paper I. We briefly recall a few relevant in- 
puts: 

- Overshooting: the convective-core baundaries are deter- 
mined with the Schwarzschild criterion. The convective core 
is extended with an overshoot pa rameter d ovel /Hp - 0.10 
(after calibration, as explained in lEkstrom et al.ll2012l) . 

- Radiative mass loss: the more massive models (M > 5 Mq) 
lose mass through li ne-driven winds, ac cording to the mass- 
loss prescription bv lde Jager et al.l(ll988l) . Only stars with an 
initial mass greater than 12 Mq lose a non-negligible amount 
of mass through radiative winds during the MS phase. For in- 
stance, our 15 Mq stellar model at Z = 0.014 loses about 
10% of its initial mass during the MS phase through ra- 
diative winds. For the red-giant p hase, the mass loss ra tes 
are obtained by the pre scription of | Reimersl(ll975l Il977l) up 
to 12 Mq, whereas the Ide Jager et al.l 0988) prescription is 
used for masses of 15 Mq or larger. 

- Mechanical mass loss: when overcritical velocities are 
reached at the equator of a star, some matter is released and 
presumably launched along Keplerian orbits in a disk around 
the star. In the following, we call such an event an overcrit- 
ical episod^3- Because of numerical difficulties in comput- 
ing of the stellar evolution when the model is at the critical 



4 It is important to properly distinguish between V/V alt and Q./Q. crit . 
Both quantities indicate how close to its breakup limit a star rotates, but 
they are not equivalent. Writing the angular velocity as SI = V cq //? e q, 
where is the rotational velocity at the equator and R ci> is the equato- 
rial radius of the star, it turns out that fl/f2 crit = (V„ /Vdt)- (Reqfi&l '^eq)- 
Because R ct] depends on the rotational velocity of the star, Sl/Sl CI it is 
equal to V e q/Krit only when SI = or in a star rotating at the critical 
limit i?cq,ciit = Req ■ Throughout this work, every time we refer to rota- 
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limit, we need to define a value for the highest fraction of 
the critical angular velocity allowed, w max , above which the 
mechanical mass loss is switched on. In our computation, we 
set w max = 0.99. The quantity of mass ejected mechanically 
from the model is governed by the quantity of angular mo- 
mentum that needs to be removed from the surface to main- 
tain the model below the limit set by w max , as explained in 
Paper I. The mass injection into the disk allows the star to de- 
crease its surface velocity below the critical value, until the 
secular evolution (and more precisely the transport of angu- 
lar momentum by meridional circulation and shears) drives 
the surface to become overcritical again. Thus, once the crit- 
ical rotation is reached, the surface velocity of the model will 
stay around the critical value until the end of the MS, i.e. un- 
til the evolutionary timescale becomes too short to allow the 
transport of angular momentum to the surface. 
- Rotation: the transport of angular momentum and chemical 
species inside a star is implemented following the prescrip- 
tion of lZahnl(ll992h for the horizontal diffusion coefficient, 
and that of lMaederl (fl997) for the shear diffusion coefficient. 

3. Properties of stars at the critical limit 

In Fig. [U we show the evolution of the co = Q/O cl -j t ratio for the 
9 Mq models at th e SMC metallicity, together with the corre- 
sponding models of Ekstrom et al. (2008, dotted lines). 

We recall that we start the computation of a mode l with solid- 
body rotation on the ZAMS, similarly to the work bv lHeger et al.l 
(2000). The underlying assumption is that pre-MS stars are fully 
convective, and that convection drives rigid-body rotation. Once 
the solid-body constraint is released, the flat internal profile of 
rotation is modified very rapidly under the action of meridional 
currents, and the surface rotation decreases strongly, as can be 
clearly seen in Fig. [TJ This means that, in the present scenario, 
the models begin their evolution with a slow-down of their sur- 
face rotational velocity, and may reach the critical limit only af- 
ter some evolution on the MS. As a result, the present way of 
initiating the computation makes it impossible to obtain very 
rapidly rotating stars at a very young age, even when impos- 
ing a high value of (0^. The validity of the assumption made is 
now being explored with models of pre-MS evolution and will 
be the subject of a future publication. The first results show that 
the internal rotation profile depends to some extent on the way 
the star mass and angular-momentum content build up through 
accretion. At the moment, though, we take the present results 
at face value, keeping in mind that we might be missing some 
young critical rotators. 

From Fig.Q] we also see that the critical velocity is reached 
at the end of the M S phase. This occured in the models by 
Ekstrom et al. (2008]) and is confirmed by the present models, 
although significant differences between these two sets of mod- 
els appear for the most rapidly rotating models. The change in 
the prescription for the diffusion coefficient explains the slight 
changes in the MS lifetimes, but has no impact on the evolu- 
tio n of the surfa c e velo city. This is clearly shown in the work 
by iMevnet et al.l d2013l) . who compared the time-averaged ve- 
locities during the MS phase obtained using different sets of 
diffusion coefficients, D,/ J(Mr and D/,. These authors found that 
for a 15 M Q model at Z = 0.002 with V;„;/V cr!t =0.5 the dif- 
ferences caused by using the present prescription and that of 
lEkstrom et all (2008) are very small, of the order of only 0.5 %. 

tion close to the critical limit, we implicitly assume that we talk about 
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Fig. 1. Evolution of the Q/Q crit ratio for 9 M© at Z = 0.002, 
with initial values of Q ini /Q cr it = 0.10, 0.30, 0.50, 0.70, 0.80, and 
0.90 (bottom to top) during the MS phase. The s olid lines show 
the pre sent models, the dotted lines the models by Ekstro m et al.l 
(120081) . 

The source of the difference in surface velocities for the most 
rapidly rotating models is the improved numerical procedure for 
computing the surface velocity. The present models maintain a 
lower surface velocity because we now account for the angu- 
lar momentum in the envelope and for the deformation of the 
star, and we also ensure the angular-momentum conservation of 
the star+wind system (including radiative and mechanical mass 
loss). As a result, the present model s evolve less easily to wards 
the critical limit, while the models bv lEkstrom et al.l d2008) over- 
estimated the surface velocity. 

The main properties of the models reaching the critical ve- 
locity at their surface during the MS phase are summarised in 
Tabled] 

From Table Q] we see that only models with <jj; n i > 0.80 - 
0.90 are able to reach the critical velocity during the MS phase. 
If we compare these limits with those given by Ekstro m et al.l 
(120081) . we see that for the metallicity and the mass range in 
common between the two works, namely for the 3 and 9 Mq at 
Z = 0.002, the prese nt models require a m inimum = 0.95 
, while the models by Ekstro met al.1 J2008) require cjm = 0.90 
and 0.80, respectively. This reflects the fact mentioned above 
that the present models approach the critical limit less easily. We 
also see that the lower-mass models reach the critical rotation 
much more easily than the more massive ones and can spend up 
to 30% of the total MS time at the critical velocity. This trend 
was also obtained in Ekstrom et al . (2008) at Z = 0.002. This 
comes from the longer time available in lower-mass models, for 
meridional circulation and shears to bring angular momentum 
at the surface. Moreover, in the higher initial-mass models, mass 
loss by stellar winds may remove some angular momentum from 
the surface and thus keep the star away from the critical limit 
(although this effect is quite weak at low metallicity). 

3.1. Surface velocities and abundances for very rapid 
rotators 

For the 9 Mq models at Z©, the minimum average velocity 
needed to reach the critical velocity is around 330 km s _1 . At 
Zsmc this minimum value is higher by about 20% (~ 390 km 
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Table 1. Data for the models that reach the critical velocity during the MS. 
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Notes. Initial metallicity, mass and rotation rate (columns 1 to 3), time spent at the critical limit in Myr (column 4) and in units of the MS 
lifetime (column 5), central hydrogen mass fraction when the critical limit is reached for the first time (column 6), total mass lost during the whole 
computed evolution (column 7), mass lost through stellar winds (column 8), mass lost mechanically at the critical limit (column 9), mechanical 
mass-loss rate averaged over the critical-rotation period (column 10), mean equatorial velocity during the MS (column 11), and enhancement 
factor of the N/C mass-fraction ratio between the ZAMS and the end of the MS (column 12). 



s~'). For the 1.7 Mq models, the minimum average velocity 
amounts to 240 and 270 km s~' at Zq and Zsmc- 

At the end of the MS phase, quite large enhancement factors 
of the N/C abundance ratios are obtained at the surface of the 
initially-rapid rotators. While at Zq it amounts to a factor of 2 
for the 1 .7 Mq model, it amounts to a factor of 20 or more for the 
9 Mq model. The corresponding enhancement factors at Zsmc 
are between 20-40 for the 1 .7 Mq model and around 500 for the 
9 Mq model. 



It is interesting to note that when these very rapid rotators 
evolve into the red-giant or red-supergiant stage, they keep some 
traces of their former very high rotation rate (see Tables 2, 3, 
and 4 in Paper I). This is reflected in two observable quantities: 
1) their rotation rate is higher, typically the model of 2.5 Mq at 
Z = 0.014 with a>i„i = 0.10 has a rotation period P mt = 950 days 
at the end of the core He-burning phase, while the corresponding 
model with a) m i = 0.95 has P mt < 150 days at the same evolu- 
tionary stage, i.e. it rotates at least six times faster; 2) the N/C 
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ratio is enhanced in formerly rapid rotators, typically the 2.5 Mq 
originating from the <x>; n i — 0.10 model has a N/C ratio of 1.4 (in 
mass fraction) at the end of the core He-burning phase, while the 
corresponding model with io m \ = 0.95 has N/C = 4.8, i.e. more 
than three times higher. These differences are stronger at lower 
metallicities. Typically, the rotation for the 2.5 Mq at Z = 0.002 
at the end of the core helium-burning phase is enhanced by a 
factor of 7 between the io m \ = 0.10 and 0.95 models. The N/C 
ratio is more than 12 times higher in the initially most rapidly 
rotating model than in the slower one. 

From the above results, we see that it might be very interest- 
ing to study red giants in clusters that show a large number of Be 
stars on the MS band. From the present models, one expects that 
the red giants that evolved from the Be-star population present 
surface velocities and N/C ratios higher than normal. 

3.2. Expected number of rapidly rotating B-type stars in 
clusters as a function of age 

The fraction of stars with surface velocities above a given limit 
varies in clusters of various ages, and it can be estimated from 
the computed grids. This is shown in Figs. [2] and [3] for solar 
and SMC metallicity. The populations were computed as fol- 
lows. First, the mass range (1 .7 Mq-15 Mq) and initial rotation- 
rate range (0 < u> m i < 1) were discretised in several intervals 
(1000 mass intervals and 100 velocity intervals), dividing the 
plane Mi n i-w;„i into 100 000 cells. For each cell, we extracted 
the expected number of stars (normalised to an initial popu lation 
of 1), according to an initial mass function (Sar peterll 19551) . and 
to an initial rotation-rate distribution dHuang et al.ll2010l) . Each 
cell thus defines initial conditions in Z, M, and w. We built an 
evolutionary track for these initial conditions interpolating be- 
tween our computed tracks. At each time step, we checked the 
spectral type of the star in each of the cells, as well as its current 
rotation rate. If this rotation rate was higher than a given con- 
stant (0.70, 0.90, or 0.98), we added the number of stars in the 
cell to the total number of stars of a given spectral type that ro- 
tates more rapidly than the given constant. The population was 
normalised with respect to the total number of B stars (includ- 
ing Be stars) at the same tim We assumed that the distribution 
of initial velocities given bv lHuang et alJ (J2010) corresponds to 
the distribution of velocities on the ZAMS, i.e. before the slow- 
ing down due to the activity of the meridional currents described 
above in Sect. [3] We emphasise here once more that an appro- 
priate account of the way the mass and angular-moment content 
of the star builds up during the accretion phase may change this 
picture. This will be studied in detail in a forthcomming paper 
(Haemmerle et al. in preparation). 

Under the previous hypotheses we find that at Zq, the frac- 
tion of stars with u> > 0.98 (very rapid rotators hereafter, blue 
solid line in Fig. [2]) becomes greater than 1% after 100 Myr, 
when the initially rapidly rotating stars with masses between 4 
and 5 Mq approach the critical limit. The fraction of very rapid 
rotators then increases with the age of the cluster until a maxi- 
mum at around 4% is reached at an age of about 300 Myr. 

The curve for u> > 0.90 (green dashed) presents a local min- 
imum around 20 Myr. This reflects the fact that in the first part 
of the evolution of our lower initial mass model, the o> ratio de- 
creases. Therefore this first behaviour is quite dependant on the 
initial conditions chosen. Then the curve follows a similar trend 
as the one for u> > 0.98, but shifted at higher values. In con- 

5 We consider a star to be a B-type star if 10, 000 [K] < T cff < 
30, 000 [K]. 




time [yr] 

Fig. 2. Evolution of the fraction of stars that rotate faster than a 
given limit in clust ers of various ages at Z = 0.014, taking the 
iHuang et al.l d2010h distribution of velocities as the initial veloc- 
ity distribution. The observational points represent the observed 
Be/(B+Be) ratio dMcSwain & Giesll2005l) . The squared symbols 
(solid error bars) are for confirmed Be stars, and the triangle 
symbols (dashed error bars) are for all Be-type stars. The ticks 
at the top of the plot show the typical lifetime of stars with vari- 
ous initial masses. 



trast, the line corresponding to w > 0.70 (red dotted) is almost 
constant up to about 400 Myr. 

Comparing the values obtained at Z = 0.014 (Fig. |2]l and 
0.002 (Fig. [3), a few interesting remarks can be made: 

- for most of the age interval between 30 and 400 Myr, the 
fraction of stars with o> > 0.98 at a given age is greater at 
Zq than at Zsmc- This occurs because 4 and 5 Mo mod- 
els at Zsmc do not attain values of to close to 1 (See also 
Table 1): even though the rotational velocities do not dif- 
fer significantly between SMC and solar metallicities for the 
mentioned masses, these models do not approach the critical 
limit because the stars are more compact and their critical 
velocities are higher at SMC metallicities. 

- Stars with u> > 0.98 are found in older clusters at Zsmc than 
at Zq. This is because at SMC metallicity, at ages older than 
1 Gyr only the lower mass rapidly rotating B-type models 
are still in their main sequence. At this age and metallicity, 
all other models with masses higher than 1.7 Mq and those 
of the same mass but with medium or low rotation, will have 
ended their core-hydrogen burning. The enhancement in the 
main-sequence lifetime for 1.7 Mo rapidly rotating models 
occurs because their central temperature attains values for 
which the pp-chain is almost entirely responsible for the en- 
ergy production during a long part of the MS phase (see dis- 
cussion and Figure 5 in Paper!])- 

- For stars with u> > 0.90, the differences between Zq and 
Zsmc ar e le ss marked, except that at the Zsmc some older 
clusters show a significant fraction of u > 0.90 stars and that 
at very high ages (above the billion years), a great fraction 
of the B-type stars may be rapid rotators. This is expected 
because only the most rapidly rotating models are still B- 



5 



Granada et al.: Populations of rotating stars 



l — i i i i I 1 1 1 — r~ 

- 1 I I'll I 

15 M„ 7 M„ 5 «„ 4 «. 3 W„ .. .A Ma 




10 8 

time [yr] 

Fig. 3. Same as Fi g. 13 for Z _ 0.0 2. Ob served Be/(B+Be) ra- 
tio extracted from Mart avan et al.l (1201 Ob . following the same 
procedure as in iMcSwain & Giesl (120051) . We consider the fol- 
lowing age beams: 7 < log(f) < 7.5, 7.5 < log(f) < 8.0, 
8.0 < log(f) < 8.5, and log(f) > 8.5. The error bars for the 
fraction of Be-type stars are Poisson statistic error bars. 



type stars on the MS band at these ages. They are some kind 
of blue stragglers^. 

Implications of these results concerning the question of the 
Be phenomenon are discussed in Sect. [6] 
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Fig. 4. Evolution of the mechanical mass-loss rate as a function 
of time for the 7 Mq stellar model at solar metallicity with a time 
step of about nine years (case n = 105, see text). Each point rep- 
resents the mass loss rate in one time step. The quiescent phase 
between two mass-ejection episodes lasts for about 40 years. 

is shown in Fig. [4] during a period spanning about 1000 years 
for our 7 Mq stellar model at solar metallicity. We see a very 
regular pattern where the quiescent phase between two mass- 
ejection episodes lasts for about 40 years, and the mass ejections 
are of the order of 10~ 8 Mq yr . For this model, the radiative 
winds are negligible and thus the only mass-loss mechanism tak- 
ing place is mechanical mass loss, as described above. 

In the following, we explore the dependency on the results 
on the time step used for the calculations and on the parameter 



4. Mechanical mass loss at the critical limit 

In this section, we explain how the mechanical mass losses (see 
columns 9 and 10 in Table Q]) were computed. We want to stress 
that in the frame of the hypotheses described below, we obtain an 
outbursting behaviour for the mass loss. However, this outburst- 
ing behaviour results from our numerical procedure and does not 
come from ab initio physical principles. In this section we ex- 
plain in detail how the mechanical mass losses were computed 
and study how the amount of mass lost depends on the time step 
and on w max . 

The general principle is the following: owing to internal 
transport of angular momentum and to the stellar evolution it- 
self, it may occur, depending on the initial rotation, that the stel- 
lar surface reaches (and even surpasses) the critical velocity (or 
rather w max in our computation, see Sect. [2} during the evolution. 
In this case, the excess angular momentum is removed from the 
outer layers in the form of an additional equatorial mass loss on 
a dynamical timescale, i.e. within one time step. This brings the 
surface velocity of the star to below 0.99 £2 c rit- Ideally, the stel- 
lar surface should rotate at 0.99 £2 cr j t , but it does not, because 
of discretisation. The transport of angular momentum by merid- 
ional currents and shears will then accelerate the surface until the 
limit is crossed again and some more mass is shed. For illustra- 
tion purposes, the evolution of the mass loss obtained this way 

6 At Z = 0.002 we slightly underestimate the number of B-type stars 
since some of them would come from stars with a mass lower than 1.7 
Mq, and thus we overestimate the Be/(B+Be) ratio. 



4. 1 . Discussion of the choice of the time step 

First, we briefly examine how the choice of the time steps modi- 
fies the amplitude of the mechanical mass-loss rates. Depending 
on 1) the duration of the time steps, 2) the equatorial velocity of 
the star just before becoming overcritical, and 3) the efficiency 
of the angular-momentum transport, the surface velocity during 
the next time step may surpass the critical limit to a lesser or 
greater extent. In some cases, very weak mechanical mass-loss 
episodes will occur (typically when the amplitude of the evo- 
lution of V eq beyond the critical limit will just graze the limit), 
while in some other cases, very strong mass-loss episodes will 
occur, when V eq strongly surpasses the limit. Therefore our nu- 
merical method leads to some scatter in the mass-loss rates. In 
Fig.|4] the interval shown is too short for this scatter to appear. 

To study the impact of the choice of the time step on the 
outputs of our models, we computed different 7 Mq models at 
Z = 0.014 with w ini = 0.95. The first column of Table in- 
dicates the value of n, characterising the different time steps 
At = A? s tandard/i) where A? stan dard = 942 yr. The last two columns 
show the total mass lost and the mean mass-loss rate during 
the considered period, which in the present case corresponds to 
a period of 250 000 yr beginning when the age of the star is 
5.21 • 10 7 yr. When the time step becomes lower than about 100 
years, these two quantities stabilise around 7.3 ■ 10~ 4 Mq for 
the total mass lost mechanically during that period and towards 
2.9 - 3.0 ■ 10~ 9 Mq per year for the time-averaged mass-loss 
rates. This shows that below a given value for the time step, the 
results are no longer varying. We also see that taking a time step 
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greater by an order of magnitude does not change the results by 
an order of magnitude: the difference amounts at most to a factor 
of two in the range of the time steps studied here. 



Table 2. Time-step dependence of the total mass lost mechani- 
cally and of the mean mechanical mass-loss rates (in the 250' 000 
yr interval). 




10- 
M [M e yr" 1 ] 

Fig. 5. Fraction of time spent with a given mass-loss rate for the 
7 Mq models at Z = 0.014 with various time steps. 



Figure [5] shows the fraction of time spent near the critical 
limit with zero and non-zero mass-loss rates for our 7 Mq model 
at Z = 0.014. Depending on the time step used, the fraction of 
the time spent away from the critical limit, i.e. when the star 
is rotating fast but does not lose mass through mechanical mass 
loss, amounts between 50 and 80% of the critical period. Shorter 
time steps tend to increase the duration of the null mechanical 
mass-loss period and to increase the mass-loss rates in the period 
when mass ejections occur. The differences become small when 
time steps shorter than 100 yr are used. 

As a conclusion, we see that the scatter brought to the results 
by a change in the time step is not dramatic. The adoption of the 
longest time steps does not change the results by more than a 
factor of two with respect those obtained with the shorter time 
steps, and also allows for a reasonable computational time for 
these extended grids. 

To support these numerical results, we compared them with 
a simple analytical estimation of the mechanical mass loss rate. 
Asuming that the star rotates as a solid body, that the rotational 
velocity change is very small when the star rotates close to the 
critical limit, and also that the time variation of the stellar an- 
gular momentum of the star is due to the loss of matter through 
the equator, we obtain simply that M = l/Rg„. This shows that 
the mechanical mass loss rate when the star rotates at the crit- 
ical limit depends basically on the evolution of the momentum 
of inertia mass and the equatorial radius, which both depend on 
the stellar mass. With this approximation, within the same in- 
terval of time as described above, we obtain a mass loss rate of 
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1.5 x 10 9 for a 7 Mq model at Z=0.014, which nicely agrees 
with the numerical results. 



4.2. Choice of the value of u> max at which mechanical mass 
loss occurs 

As emphasised above, because of the numerical difficulties of 
computing the stellar evolution when the star is at the critical 
limit, we need to define a value for the highest fraction of the 
critical velocity allowed, which we call w max (see Sect. |2]i above 
which the mechanical mass loss sets in. 

To check to what extent the choice of u> m!ix affects the evo- 
lution of the star, we computed 7 Mq models at solar metallic- 
ity with win; = 0.95, setting w max equal to 0.99, 0.95, and 0.90 
just after the readjustment phase at the very beginning of the 
MS. The results are presented in Table|3] The first column gives 
w mas , the second and third one the MS lifetime and the fraction 
of the MS lifetime spent near the critical limit; the total mass lost 
and the time-averaged mass-loss rate over the critical period are 
given in columns 4 and 5. 

We see that the lower w max , the longer are the MS lifetimes, 
and the higher it is the ratio T cr i t /TMS- This is a consequence of 
the surface braking by the mechanical mass loss, which occurs 
earlier for the models with the lower values of w max , producing 
a stronger internal mixing in the star. This mixing brings fresh 
hydrogen to the core, increasing its size, slightly increasing its 
luminosity, and also the MS lifetime. The total mass lost and 
the mechanical mass-loss rates are also higher for lower w max . 
This is expected, since a lower value of w max implies an ear- 
lier entrance into the critical regime, more mass removed when 
the limit is crossed, and a higher amount of total mass lost. 
Typically, passing from w max = 0.99 to 0.95 already increases 
the total mass lost by more than a factor of two. Because both 
the total mass lost and the duration of the critical period increase 
in a similar way when io^ decreases, the averaged mass loss rate 
remains nearly constant. 

Figure [6] shows the HR diagram for the 7 Mq models com- 
puted with the different values of w max , and as a comparison the 
evolutionary tracks corresponding to the non-rotating case and 
the evolutionary tracks corresponding to <w max = 0.99, for 5, 7 
and 9 Mq. The differences in log(L) and log(r e ff) have different 
causes. As mentioned above, the luminosity shift is due to the 
larger core of the model with the lower <w max . The T e s shift has a 
geometrical origin. For a given luminosity, the surface of a star 
rotating with o> = 0.99 is 1.16 times larger than the surface of 
a star with oj = .90, because of the deformation of the shape 
of the surface rsee l Geomv et al. lEoTTh . Because the mean T c q is 
deduced according to the Stefan-Boltzmann law L 



zr 4 ff (with 
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Fig. 6. HR diagram for 7 Mq models at solar metallicity com- 
puted with different values of co max (blue: 0.99, green: 0.95, and 
red: 0.90). As a comparison 5 and 9 Mq evolutionary tracks are 
also plotted. 

Table 3. Properties of the 1 Mq stellar models at Z = 0.014 
computed with various values of <y max - 
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2 the stellar surface), this produces a shift of 0.016 dex towards 
the red in the HR diagram. 

In the rest of the paper, we discuss the results obtained with 
<^crit = 0.99. To conclude this section we stress two points: 

1. Reality may be quite different from the situation described 
above. It might occur, for instance, that mechanisms other 
than critical rotation are involved in producing mass ejection. 
These additional mechanisms may allow the loss of matter 
well before the surface rotational velocity becomes critical. 

2. While the present approach is of course probably too 
schematic, it can at least provide a reasonable estimate of 
the mass lost mechanically when the following hypotheses 
are made: mechanical mass loss only occurs when the veloc- 
ity is higher than a fixed limit and the quantity of mass lost 
is entirely determined by the necessity of the star to shed 
overcritically rotating material. 



4.3. Radiative and mechanical mass loss 

In Table [TJ column 10, we show the mass-loss rates averaged 
over the period during which the star evolves near the critical 
limit for stars of different stellar masses and metallicities. 

The total mass and angular momentum lost mechanically 
are quite small with respect to the total mass and angular- 
momentum content. Typically, for the 15 Mq model at Z = 



0.002, the total mass lost mechanically corresponds to 0.005% 
of the initial mass. This means that this mass-loss mechanism 
has only very little influence on the evolution of the star. 

For a given metallicity, the mean mechanical mass-loss rates 
are higher for higher stellar masses, As mentioned in subsection 
4. 1, this dependence can also be deduced from simple analytical 
calculations: the mass-loss rates depend on the evolution of the 
stellar momentum of inertia and on the equatorial radius of the 
stars, and both quantities depend on the initial stellar mass. For 
a given mass, we find that the mechanical mass-loss rates are in 
general higher at lower Z. In the high mass star range, this can be 
explained easily: indeed, in the high-mass range, the stronger ra- 
diative winds of higher Z remove mass and angular momentum 
at the surface, keeping the star more easily away from the criti- 
cal limit and thus reducing the need for a mechanical mass loss. 
In the low-mass range, the behaviour is more complex, which is 
the result of counteracting effects. First, stars are more compact 
at low Z, so the extent of the region over which angular momen- 
tum must be transported is smaller; second, this compactnes s 
leads to weaker meridional currents (Maeder & Mevnet 20011) . 
which reduces the transport efficiency. The present results seem 
to imply that in the metallicity range presented here, the first ef- 
fect overcomes the second. From the above, we conclude that 
our mechanical mass loss rates increase with mass, as is the case 
of mass loss rates via radiatively line-driven winds, but mechan- 
ical mass loss is higher at lower Z, in contrast to what happens 
for stellar winds. 

An interesting point to be highlighted here is the impact of 
the mass loss on the angular-momentum content of the stars. The 
total initial angular momentum of models that rotate with differ- 
ent w; n j remains within an order of magnitude for a given stellar 
mass. Along the main-sequence evolution, the angular momen- 
tum of the star will be modified by the effect of stellar winds 
(for stars with masses larger than 7 Mq) and by the mechanical 
angular-momentum loss, for stars that reach the critical limit. 
The stellar winds dominate the high-mass range, driving the loss 
of about 30% of the initial angular momentum at Z = 0.002. The 
mechanical mass loss reduces the angular momentum at most up 
to 10%, but is very effective in injecting angular momentum into 
the circumstellar environment, probably forming a circumstellar 
disk. 



4.4. Polar and equatorial mass fluxes 

If mass loss occurs only via radiatively driven winds, the equa- 
tor ial to polar mass-flux ratio depen ds only on the value of 
a) (|Maeder 2002t iGeorgy et al.ll201 lb and basically reflects the 
temperature contrast between the pole and equator described by 
the von Zeipel's theorem. For a star rotating at u> = 0.99 the 
equatorial to polar mass-flux ratio is therefore pr edicted to be 
at mo st 0.15, if there is no mechanical mass loss (Georg y" ;t al.l 
120111) . We see in Table g] that for the models with 7 and 9 Mq, 
the mean mass flux is higher at lower metallicities. This might 
be surprising at first sight since a stronger mass loss by radiative 
winds is expected at high metallicity because of the metallicity 
dependence of the line-driven winds (which scale as ~ Z" 5 ). 
However, polar winds occur only when the surface velocity is 
close to the critical limit and the evolutionary stage when this oc- 
curs depends on the metallicity, and at low metallicity the critical 
limit is reached closer to the end of the MS, where the radiative 
winds are stronger due to the higher luminosity of the star. 

When mechanical mass loss occurs, the mass flux at the 
equator consequently becomes much greater than the mass flux 
at the pole (enhancement factors of up to 60) in most cases. 
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The only exceptions occur for the highest initial masses in each 
metallicity. In that case, the radiative winds are stronger and thus 
remain of the same order of magnitude (or even stronger in one 
case) than the equatorial mechanical mass flux. 



5. Decretion disks 

5.1. Estimation of the physical properties 



IStruvd (Il93ll) first proposed the existence of a decretion disk 
formed from an equatorial-mass ejection episode by a rapidly 
rotating star. Since then, the existence o f Keplerian rotating disks 
has been incr easingly supported (e.g. iPorter & Riviniusll2003t 
ICarciofill20Tll) . 

In the previous section, we showed that the surface veloc- 
ity remains at around the critical limit from the first moment it 
reaches it until the end of the MS phase, and that mass ejec- 
tions in the frame of the hypotheses made here occur as inter- 
mittent outbursts. The ejected mass launched into an equatorial 
disk around the star explains the observed Balmer emission lines 
of Be stars. 

We study here the physical characteristics of a station- 
ary disk that can be deduced from the developments made by 
iKrticka et al.l d201 ll) . who solved the hydrodynami c equations in 
cylind rical coordinates to obtain the disk structure. Krtick a~et al.l 
d201 ll) provided parametrised expressions for the mass-loss rate 
as a function of the angular-momentum-loss rate and for the 
outer radius of the disk. We used these simple parametrised 
expressions to study how the disk properties change along the 
critical-rotation phase. 

There are three free parameters in this disk model: two of 
them, p and To, arise from the assumption that the temperature 
in the disk varyies as a power-law in radius, T - To(R eq /r) p , 
where /? eq is the equatorial radius of the star. The third f ree pa - 
rameter is the viscosity term a of Shakur a~& Sunvaevl (Il973l) . 
which intervenes in the momentum equation (tangential compo- 
nent). This viscosity couples the inner parts of the disk to the 
central star and may be generated b y the presence of small mag- 
netic fields dBalbus & Hawlevlll991l) . 

The structure of the disk and our stellar model are connected 
through three variables: 1) X, which is the quantity of angular 
momentum that the star must lose per unit of time to maintain 
the surface velocity at the critical limit (at w max = 0.99); 2) R eq , 
the equatorial radius of the star; and 3) VkCReq)> the Keplerian 
velocity at the equator of the star (see more details below). 

Part of the excess angular momentum that the star has to lose 
is lost through radiative stellar winds. To calculate the angular- 
momentum-loss rate by this process, -Cwind, we inserted the ra- 
diative mass-loss rates M given by our evolutionary models and 
assumed that this mass is lost isotropically at the mean stellar 

radius given by R mesin = ^L/AmrT^, which gives £ wind = 

2/3MilR^ ean . By doing so, we neglected the wind anisotropies 
and thus slightly overestimate the loss of angular momentum due 
to this process. However, except for the 15 Mq stellar model 
at Z = 0.002, the stellar winds remove a negligible fraction of 
the mass and of the angular momentum even when the effects 
of the stellar winds are overestimated. The remaining angular- 
momentum-loss rate, £, - -Cwind, is transferred to the disk (£, 
hereafter for clarity). 

In a stationary situation, the disk will lose angular momen- 
tum at the same rate. Therefore, knowing how and where the 
disk loses this angular momentum will set up the disk structure 



corresponding to the considered physical situation. The disk can 
lose angular momentum through two different channels: 

1. through disk winds (_£dw)- The value for t he maximum 
. £dw i s obtained from expression (26) in Krticka et al. 
(1201 ll) . which corresponds to the highest disk-wind angular- 
momentum-loss rateQ In all our models the highest _£dw is 
negligible with respect to £. (compare column 5 with col- 
umn 3 in Table 0, therefore we neglect the effect of a disk 
wind hereafter; 

2. through disk mass loss (-Cdisk)- Since this is the only active 
process in most cases, one has X = -Cdisk- 

The quantity £. is given by our evol utionary calculations , 
while the value of -Cdisk is obtained by Krtic ka et alj d201 lb . 
These authors showed that at a radius R a j t the radial velocity 
of the matter in the disk becomes equal to the local sound veloc- 
ity, a - -JkT/iimu, where fi is the mean molecular weight and 
»ih the mass of a hydrogen atom. Matter flowing beyond R cr n 
removes angular momentum from the star-disk system. 

Accounting for the fact that at large distances the disk is not 
rotating at the Keplerian velocity, 

-Cdisk = -^diskVK(^crit)^crit, 

where Md; s k is the mass lost by the disk, VK{R C nt), the Keplerian 
velocity at R CI \ t , i.e. Vic^nit) = y/GM/R a i t . The equality be- 
tween X and Xdisk allows us to determine Mdisk, provided R CT \ t is 

kn own. 

|Krticka et al. (201 1) proposed an approximation of the crit- 
ical radius, which can be considered as the outer radius of the 
disk: 



Rout 
Req 



Rr 



VK(Req) 



10 + 4/? \ a(R eq ) 



(1) 



where a(R eq ) is the sound speed evaluated at R eq . 

At this point, we have three interesting quantities that char- 
acterise the disk: its radial extension, R a j t (defined here as its 
extension up to the point when angular momentum is no longer 
conserved), the rate of angular-momentum loss, .Cdisk, and of 
mass loss, Miisk- 

We easily obtain two more interesting quantities, the diffu- 
sion time and the mass of the disk. The characteristic timescale 
for the evolution of the disk is given by the viscous diffusion 
time, 



fdiff - 



R 1 < 

cnt 



(2) 



where v = aaH is the diffusivity of the disk, with a 
the viscosity parameter, and H the disk scale height, H = 
a(R)R^t 2 /(G M) 1 ^ 2 , where G is the gravitational constant and 
M the mass of the star. Evaluating H at R ct n and assuming an 
isothermal disk leads to 



fdiff - 



"JGMR C 



(3) 



From the diffusion timescale and mass-loss rate of the disk, 
we can estimate the mass of the decretion disk. It is obtained 



7 If this highest _£ dw were be greater than the angular momentum car- 
ried by the disk outflow, X- Xwind, the extension of the disk would be set 
by the radius R oul , given by Krticka et al. (201 1) such that -Cdw(^out) = 
-t - -Cwind- This never happens in our models, see text. 
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Table 4. Mean mass fluxes during the critical-rotation phase for stars with Quu/fi C rit = 0.95. 
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Notes. Mean polar and equatorial mass fluxes due to radiative winds (columns 3 and 4 respectively), mechanical mass flux (column 5), polar to 
equatorial mass-flux ratio due to anisotropic winds (column 6), and polar to mechanical mass-flux ratio (column 7). We only consider here models 
for which we included a radiative mass-loss prescription during the MS (M ini > 7 Mq). 



6000 




10 

t [10 7 yr] 

Fig. 7. Evolution of the radius and diffusion time of the disk as 
a function of time for critically rotating models with an initial 
mass between 4 and 9 Mq at solar metallicity. 



by multiplying the mass-loss rate of the disk by the diffusion 
timescale: 



M d: 



2l 



isk 



VGMfi c 



T=. (4) 



The mass of the viscous decretion disk does not depend on the 
outer radius R a i t but is set by the angular-momentum-input rate 
in the disk when the stellar surface reaches the critical limit. 

The disk properties are computed assuming an isothermal 
disk (p - 0) of temperature Tq - 2/3 (which are the disk 
parameters frequently used in t he literature to repres ent Be-star 
disks), and a=l as proposed bv lCarciofi et al.l (120121) . 



5.2. Results from the stellar models 

In Table |5J we give the properties of the disks that are obtained 
by applying the equations of the previous section to our stellar 



models with to^ = 0.95. The indicated quantities are obtained 
by computing their time-averaged value throughout the critical 
period to remove the scatter due to the outbursting nature of the 
process. 

We note the following points: 



- The mass-loss rates via the viscous decretion disk, (Mdisk), 
are much lower than the mass-loss rates given in Table [TJ 
A^mech, mean, corresponding to the mass lost at the equator of 
the star: in the disk, the matter is released at a farther dis- 
tance, and thus a lower mass loss is needed to remove a given 
amount of angular momentum. The mass-loss rates from the 
disk are about a fa ctor of 10 below the mass-loss rates at the 
surface of the star (Krticka et al. 2 01 ll) . 

- The masses of the disks are quite small. Values between 9.4 • 
10~ 12 and 1 .4- 10~ 7 M Q (310 6 to 0.047 times the mass of the 
Earth) are obtained. These quantities are strongly dependent 
on the disk parameters p, To, and a. 

- Under the hypotheses considered, the extension of the disk 
amounts to a few thousands of solar radii (between about 
10 and 30 Astronomical Units). Expressed in terms of stel- 
lar radii, the extension amounts to some hundreds of stellar 
radii. 

- The diffusion time is very short, only a few decades. This is 
much too short for any instability in the disk to set up and to 
form something like an asteroid belt! 

- At a given metallicity, the mass and the extension of the disk 
increase with the initial mass (see upper panels of Figs.[7Jand 
|8}. Typically at solar metallicity, the mass of the disk for the 
9 Mq model is a factor of 1600 higher than that of the 1.7 
Mq model, while the extension of the disk is increased by a 
factor of 2.2 only. Therefore, the disks are expected to be on 
average denser around more massive stars. 

- For a given initial mass, the extension of the disk increases 
when the star evolves along the MS as shown in Figs. |7]and 
[8] The same is true for the disk diffusion time. 

- For a given initial mass, in general the mass of the disk 
slightly increases when the metallicity decreases and the ex- 
tension of the disk slightly decreases. For instance, the disk 
of the 9 Mq model is 2.8 times more massive at Z = 0.002 
than at Z = 0.014, while the extension of the disk, in con- 
trast, is 16% smaller at Z = 0.002 than at Z = 0.014, i.e. 
some denser disks are expected in low-metallicity regions. 
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° 3000 




Table 6. Mean values of some disk's characteristics for the 7 Mq 
model at solar metallicity, varying Tq/T^, p, and a. 
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t [10 7 yr] 

Fig. 8. Same is in Fig. Q for an initial mass between 1.7 and 3 

Mq. 



Overall, Table [5] shows that the differences due to the metal- 
licity are small compared to those due to the initial mass. 
- The disk diffusion time does not present very strong varia- 
tions either as a function of the initial mass or as a function 
of the metallicity. 



5.3. Sensitivity to the disk parameters 

It is important to check to what extent the choice of the param- 
eters that describe the thermal structure of the disk and the vis- 
cosity parameter affects our results. We use the data of the 7 Mq 
model at solar metallicity, two different values for the parame- 
ters T /T eS (0.66 and 0.80), for p (0.0 and 0.1), and for a (0.1 
and 1.0). The results are summarised in Table [6] 

A scatter is produced in the mass, extension, diffusion time, 
and mass-loss rate of the disk. The highest-to-lowest value ratio 
varies by a factor of 33 for the mass, of 2.3 for the extension, of 
20 for the diffusion time and of 1.4 for the mass-loss rate. This 
illustrates the absolute values obtained in Table[5]are sensitive to 
the physical ingredients of the disk models. 

The mass and the diffusion time of the disks are the most sen- 
sitive quantities, and a, the viscosity, is the most critical quan- 
tity. A stronger viscosity implies smaller masses for the disks 
and shorter diffusion times. 



6. Possible link between very rapid rotation and the 
classical Be-type star phenomenon 

6.1. Is critical rotation a necessary and sufficient condition 
for Be-type stars? 

There is no doubt that Be stars are rapid rotators, i.e. with 
co > 0.70, however, it is not clear whether Be stars are ro- 
tating at the critical l imit o r not. In the sample studied by 
iLevenhagen & Leister! (|2004), eight Be-type stars with a mass 
estimate between 8 and 10 Mq show V sin i values between 160 
and 335 km s _1 , with a mean value of 230 km s _1 . Corrected for 
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0.1 
0.1 



0.1 
1.0 
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1.0 
0.1 
1.0 
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1.28 
1.28 
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1.41 
1.41 



10 -io 
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10 -io 



4.43 ■ 
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1.59 
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1.34 
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Notes. Disk parameters (columns 1 to 3), diffusion time (column 4), 
outer radius (column 5), mechanical mass-loss rate (column 6), and 
mass of the disk (column 7). The quantities are time-averaged over the 
critical-rotation phase duration. 



a random orientation of the rotational axis, this value translates 
into a mean V eq around 290 km s _I . The critical velocity for our 
9 Mq stellar model at the end of the MS phase (when the most 
rapidly rotating models reach the critical limit) is 410 km s -1 . 
Comparing these two values, we obtain a mean ratio of V/ V c rit = 
0.71 (i.e. co = 0.88) for the observed stars, which means that 
Be-type stars may not necessarily be at the critic al hmit. The 
same conclusion can be deduced from the work of iFremat et al.l 

(2005) , who found a mean ratio of co = 0.88 from the fit of 
photospheric lines of Be stars, accounti ng for gravity darken- 
ing, w hich also agrees with the result of iLevenhagen & Leister! 
(2004), who did not account for it. Be-shell stars are assumed to 
be Be stars seen equator-on. Thus, the determin ation of V sin / 
for the se objects gives direct information on V Pq . iRivinius et al.l 

(2006) showed that for the brightest shell stars of different spec- 
tral types, the mean fraction of the critical velocity is 8 1 %, which 
corresponds to a mean co higher than 0.92. Meill and et al.l(l2012l) 
found a mean value of co — 0.95 for a sample of eigh t Be stars, 
considering V sin / and V cr j t from lFremat et al.l (12 005 ) but ; mea- 
surements derived from interferometry. iDelaa et all (1201 It de- 
duced from interferometric observations of Be stars 48 Per and 
*P Per that both stars rotate at nearly critical rotation. Therefore, 
the question to which extent Be stars are really critically rotat- 
ing stars is not settled yet. Two additional factors have to be ac- 
counted for in that discussion: 1) when the surface velocity ap- 
proaches the critical limit, the method of velocity measurement 
based on the observed width of the absorpt ion lines might un- 
derestimate the values dTownsend et al.ll2004l) : 2) it may be that 
Be stars are not always critically rotating, and undergo periods 
in the course of their evolution during which they rotate below 
the critical value and periods when they are rotating critically. 

In Fig.[2]and[3] we show the fraction of Be stars observed in 
clusters wi th ages between 10 and 300 Myr at solar metallicity 
(data from McSwai n~& Gie si 120051) and betwee n 10 Myr and 1 
Gyr a t the SMC metallicity (data extracted from Mart avan et al.l 
l2010h . and compare them to the ratio given by our models. 
Considering only single stars with co > 0.98, the numbers are 
much too low to account for the observed fraction of Be-type 
stars, particularly fo r young clusters. Th is last conclusion agrees 
with the study by lEkstrom etaD (120081) . Stars with co ~ 0.80 
cannot reproduce the observed trend with age either. Then, in 
the frame our hypothesis, our present models do not provide a 
reasonable fit to the observations. 
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Table 5. Main characteristics of the decretion disks for models with Q. m \/Q. cr n = 0.95. 
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Notes. Metallicity of the star (column 1), stellar masses (column 2), total angular-momentum loss rate at the critical limit (column 3), angular- 
momentum loss rate via the stellar wind (column 4), angular-momentum loss via the disk wind considering R ou[ — > oo (column 5), angular- 
momentum loss rate via the viscous decretion disk (column 6), mass-loss rate via the decretion disk (column 7), extension of the disk expressed in 
units of Rq (column 8) and stellar radius R+ (column 9), diffusion time considering a= 1 (column 10), and total mass of the disk (column 1 1). In 
all cases an isothermal disk was considered (p=0) and a value of a=l. 
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Fig. 9. Evolution of the Be/(B+Be) ratio in clusters of various 
ages at Z = 0.014 obtained from our models assuming a mass 
(or T e g) depend ence for the Be phenomenon as proposed by 
ICranm er (2005J using (Veq/VcriOmin values). The observational 
points are the same as in Fig. [2] 



A different conclusion is obtained when one considers that 
the lowest velocity that a star has to reach to present the Be 
phenomenon may vary as a function of t he stellar mass a nd/or 
the temperature. T his was suggested by Cranmer (2005]) and 
iHuang et all d2010h . who found that more massive (or hotter) 
stars become Be-type stars at lower V/Vait values than lower- 
mass (or cooler) stars. For instance, if we consider the lowest 
rotation rate va lues (V e q/Vcrit)min in bins of T e g proposed by 
ICranm er (2005) to obtain a Be star, then the fractions of stars 
with velocities higher than these values (which appear as Be 
stars) given by our models reproduce the observed fraction of 
Be stars quite well (see F ig. 13. Our models accordingly support 
the analysis by ICranmerl (120051) . 

In addition to the variation of the fraction of Be stars with 
age, another important feature is the way this fraction varies 
as a fu nction of metallicity at a given age. It has been e stab- 
lished (Maede r et alj|1999t FvVisniewski & Bjor kmanll2006i see 
also our Fig. [2] and [3) that there is a larger fraction of Be stars 
at low metallicity. Indeed, the observed fraction of Be-type stars 
in clusters with ages between 10 and 25 Myr is higher by about 
a factor of 2 at SMC metallicity (Be/(Be+B) fraction of about 
30%) than at solar metallicity (15%). Starting from identical ini- 
tial velocity distribution, we did not obtain such an enhancement 
in our models. Solutions of this problem can come from different 
directions. For instance, it may be that the initial rotation distri- 
bution is shifted to higher velocities at low metallicity. Another 
possibility could be that since disks are denser at low metallicity 
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(for a given mass and age), the y could be more easily detected 
at low metallicity. According to lSilai et al.l (l2010). the minimum 
density for a disk to pr oduce Ha emis sion is 5 ■ 10~ 13 g cm" 3 . 
Finally, the analysis of Cranmer (2005), which sets a minimum 
velocity for becoming a Be-type star that varies as a function of 
the effective temperature was performed at solar metallicity. At 
present we cannot say whether such an analysis would hold at 
lower metallicities and whether this would help in resolving the 
above question. 

On the basis of our single-star models, we therefore conclude 
that there are much too few stars that rotate strictly at the criti- 
cal limit to account for the observed Be-type star population. If 
the Be phenomenon appears at lower rotation rate, as suggested 
by many authors, a nd has a depend enc e on the mass and/o r the 
r e ff, as proposed bv lCranmerl (120051 ) and lHuang et al.l(l2010l) . the 
models reproduce the fraction of Be-type stars much better. As 
a consequence, some other mechanism is required in addition to 
the centrifugal force to push the matter into a disk. 

This conclusion of course depends on the physics included 
in the models. Here, we did not account for any internal dy- 
namo, which would impose a solid-body rotation during the MS 
phase and in turn would generate higher surface rotations dur- 
ing the MS phase. Moreover, in addition to the mechanism ex- 
plored here, where the surface is accelerated by internal pro- 
cesses, some stars may acquire rapid rotation by tidal interac- 
tions or through mass exchange in close binaries. Furthermore, 
as already discussed previously, we started the computation on 
the ZAMS with a flat £2 profile, which implies an internal read- 
justment period at the very beginning of the evolution. As a re- 
sult, the effective Wi„i is lower than the one by which we label our 
models, and we do not produce critical rotators on the ZAMS. 



6.2. Comparison of the mass-loss rates 



IWatersI (119861) and lWaters et al.l(ll987l ) studied the density struc- 
ture and mass-loss rates of Be stars and interpreted the observed 
IR excesses in terms of a simple pole-on disk model, which al- 
lowed them to estimate upper limits for the IR (equatorial) mass- 
loss rates between 7 ■ 10~ 9 and 2 ■ 10~ 8 Mq yr _1 , and masses 
around 1 - 5 ■ 10~ 8 Mq. The thin dotted line in Fig. \W\ shows 
this u pper limit for IR mass-loss rates obtained by Waters ~et al.l 
( 1987), and the colour points in this figure show some of their 
data for which relia ble stellar-m ass determinations are avail- 
able in the literature. ISted (120031) also obtained mass-loss rates 
of some Be stars with SIMECA, a latitude-dependent radiative 
wind plus dense disk model that produces both spectroscopic 
and interferometric synthetic data, plotted in crosses in Fig. [TUl 
ICarciofi et al.l (1201 2|) obtained similar values for 28 CMa by 
modell ing the decline rate of th e V-band ex cess. Since the val- 
ues bv lWaters et al.1 (1 19871) and lSted d2003) are strongly model 
dependent, they cannot be taken as strong constraints of the mod- 
els. However, it is interesting to compare them with the mass loss 
rates we obtained in the present work. In Fig. [I0]we also plot 
the mean equatorial mass-loss rates obtained from our critically- 
rotating models as described in Section|4]and presented in Table 
Q]as a solid black line for stars with Q/Q crit =0.95 on the ZAMS. 
In this plot we also represented the averaged mass-loss rates 
for overcritical episodes (black dashed line). As mentioned in 
Section[5] the disk mass-loss rates are lower than these values. 

We see that our ma ss-loss rate estimate s run below the de- 
terminations made by Waters et al. (1987). The origin of this 
shift between the two sets of values is difficult to assess, since 
both may be model dependent. However, it is remarkable that 
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Fig. 10. Mass loss - luminosity diagram . The colour points cor- 
respond to data by IWaters et al.l (1 1987b . while the black crosses 
correspond to data from IsteeTfeOOSb and the black circle to 
the data from ICarciofi et al.1 (1201 2l) for 28 CMa. The thin dot- 
ted curvejraces t he up per limit for IR mass-loss rates obtained 
by IWaters et al.l (119871) . The thick black curve corresponds to 
the mean mechanical mass-loss rates along the critical rotation 
phase, obtained as described in Section [4] and shown in Table [1] 
while the dashed black line corresponds to the mean value dur- 
ing the time in which the star is effectively losing mass. The grey 
region corresponds to the range of instantaneous mass-loss rates. 



the models predict an increase of the me chanical loss rate w ith 
the luminosity similar to that obtained bv lWaters et al.l (119871) . 

In this context, it is interesting to note that per se, an over- 
critical velocity is not sufficient to drive a mass loss. When it 
becomes overcritical, matter will be launched into a Keplerian 
orbit and therefore will remain gravitationally locked to the star 
(Jupiter is orbiting around the Sun with a critical velocity !). This 
matter may even fall back on the star, depending on the forces 
acting on it. Some additional force (or forces?) have to come 
into play and provide the momentum needed to drive an effec- 
tive mass loss from the star and/or an outwards-expanding disk. 
Here we did not account for any additional forces. If they were 
present, they would probably increase the mass "mechanically" 
lost by a fast-rotating star, and therefore our present estimate 
must be considered as a lower limit. 

Looking at Fig. [TO] we can see that the upper values ob- 
tained in our models for the instantaneous mass loss rates fall 
in the range of the observed values. However, while the time- 
averaged mass loss rates are meaningful quantities in the frame 
of our hypothesis, this is not the case of the instantaneous values, 
which are quite model dependent. Therefore, we cannot deduce 
any firm conclusion from this coincidence. 

An interesting point is that stars rotating at or near the critical 
velocity may or may not have an equatorial disk, but all these 
stars should present some degrees of anisotropy in the radiative 
winds. As mentioned in Sect. 14.41 the radiative mass flux ratio 
between the equator and the pole reaches 0.15 at most. When 
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mechanical mass loss occurs at the equator, the ratio increases 
to 60. Some observations have been performed on a Ara (a 9 - 
10 Mq star estimated to rotate at V eq /V CT n = 0.91, i.e. u = 0.98- 
0.99) and a Eri (Achernar , a 6 - 8 Mp, st a r estim ated to rotate 
at a) — 0.992). For a Ara. IChesneau et all (f2005) found a polar 
flux of 1.7 • 10~ 9 Mq yr _1 sr _1 , with a ratio between the polar 
and equatorial flux <t> en /O nn i _ 30. These numb ers have been 
revised in a second paper (Meill and et al.l 120071) to O po i = 7 ■ 
10" 9 M q yr 1 sr _1 and O eq /<D po i = 0.03. For gEri. lKanaan etalJ 
(120081) found a polar flux of 4.0 • 10~ 12 Mq yr 1 sr" 1 and a ratio 
<t> eq /(t>p i = 0.01 - 1. The observational range is wide, and for 
the time being, it seems to be premature to use these kind of 
observations to put constraints on the models. 

6.3. Comparison of the disk mass, extension, and diffusion 
time 

In our model, we obtained that the masses of disks in crease sig- 
nifica ntly with the stellar mass. Interestingly enough, Silai et al. 
(2010) modelled Hor emission line profiles of Be stars with their 
non-LTE radiative transfer code BEDISK by fitting disk parame- 
ters: their Table 3 shows that disks around B0 type stars are more 
massive than disks around B8 type stars. 

Our models predict extensions for the disk's outer radius of 
some thousands of solar radii. Radio studies provide informa- 
tion on the outer r egions of disks surro unding Be stars. Using the 
Very Large Array. iTavlor et al.l (1l990b carried out a radio survey 
of 21 Be stars selected for their excess IR emission detected by 
IRAS. They found that the minimum radius of Be-star disks are 
of several hundre ds to thousands o f solar radii. Identical results 
were obtained bv lClark et all (119981) . who observed a sample of 
13 Be stars of the southern sky with the Australian Telescope 
Compact Array. The values of a few thousand solar radii we ob- 
tain with the different stellar masses and different disk param- 
eters seem to be compatible with the observed quantities. The 
extension of the disk may be dependent on the wavelength used 
to investigate it. In Ha, the extensions are much shorter, a few 
tens of solar radii, since only the portion which is nearby the 
star and therefore hot and dense enough would be observable, 
while at longer wavelengths, cooler and shallower parts of the 
disk can be seen. Our predictions are valuable for isolated stars. 
Disk truncation can happen if a close companion star is present. 

The diffusion timescales we obtain with our models can be 
compared to long-term variability cycles observed in Be stars. 
Some of the observed long-term variations in Be stars, consist- 
ing of transitions betwe en normal B and Be phases, or bright- 
ness variations (see e.g. Porter & Rivinius 2003: ICarciofife Ol 1. 
and references therein), can be understood in terms of the forma- 
tion and dissipation of a circumstellar disk. They typically last 
from several months or years to decades. The values of a few of 
decades we obtain are consistent with these timescales, though 
different choices of disk parameters can affect this quantity. 

7. Conclusions 

With an improved version of the Geneva stellar-evolution code, 
we have studied the characteristics of initially very rapidly- 
rotating stars as well as some properties of the equatorial disks 
that are produced when the surface velocity becomes overcriti- 
cal. 

The main conclusions are: 

- The angular-momentum and mass losses (between 4 • 10 -13 
and 4 ■ 10~ 9 Mq yr _1 , depending on the initial mass and 



metallicity) when stars are at the critical limit are very mod- 
est and have no strong impact on the way the star evolves. 

- Nevertheless, signatures of rapid initial rotation on the 
ZAMS remain present in more evolved stages such as the 
red-giant or supergiant stages, for instance the higher val- 
ues for rotation and surface nitrogen enrichment compared 
to stars that evolved from slowly rotating progenitors. 

- For stars that reach the critical limit, the time-averaged mass- 
loss rates we obtain are one order of magnitude lower than 
those currently estimated from observations of Be-type stars. 
However, among the range of instantaneous mass-loss rates 
we obtained, the highest values agree relatively well with the 
estimates derived from observations. 

- When the star is losing mass at the critical velocity, the mass 
of the disk produced is estimated to range between 9.4- 10~ 12 
and 1.4 ■ lO^M© (from 3 ■ 10" 6 to 47 • 10" 2 M e ), its ex- 
tension between 2000 and 6500 Rq, and its typical evo- 
lution timescale between 10 and 30 yr. Some hints were 
given on how these values vary as a function of the initial 
mass, age, and metallicity. With our simple assumptions, 
even though the mass-loss rates are somewhat lower than 
those derived from observations, the masses, extensions, and 
diffusion times of the disk obtained here are consistent with 
the current estimates of these quantities for observed disks 
around Be-type stars. 

- If we assume that there exists a r e (j dependence for the min- 
imum a) for the Be phenomenon to occur (as suggested in 
the literature), the models agree reasonably well with the ob- 
served values of Be/(B+Be). 

Our results were obtained in the frame of the shellular rota- 
tion model. If an internal magnetic field would more strongly 
couple the core to the envelope, and provided no substantial 
magnetic braking occurs at the surface, the star may reach (at 
an earlier evolutionary stages) the critical limit more easily. In a 
forthcoming paper we will explore the consequences of a strong 
interior coupling. 
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